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Abstract. It is said that the finding of the maximum oscillation in neutrino oscillation by Super- 
Kamiokande is one of the major achievements of the SK. In present paper, we examine the assumption 
made by Super-Kamiokande Collaboration that the direction of the incident neutrino is approximately 
the same as that of the produced lepton, which is the cornerstone in their L/E analysis and we find this 
approximation does not hold even approximately. In the Part 2 of the subsequent paper, we apply the 
results from Figures 12, 13 and 14 to L/E analysis and conclude that one cannot obtain the maximum 
oscillation in L/E analysis which shows strongly the oscillation pattern from the neutrino oscillation. 

PACS. 1 3.15.+g, 14.60.-z 
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^] ', 1 Introduction as, say, Sub-GeV e-like, Multi-GeV e-like, Sub-GeV /i-like, 

t~-! • Multi-GcV //-like, Multi-ring Sub-GeV /x-like, Multi-ring 

O \ According to the results obtained from the Super-Kamio- Multi-GeV //-like, PC, Upward Stopping Muon Events and 

kande Experiments on atmospheric neutrinos, it is said Upward Through Going Muon Events. In other words, all 

T! I that oscillation phenomena have been found between two possible interactions by neutrinos, such as, elastic and 

neutrinos, t^and v T . Published reports on the confirma- quasielastic scattering, single- meson production and deep 

tion to the oscillation between the neutrinos, z^and u T , scattering are considered in their analyses. Furthermore, 

and the history foregoing these experiments will be criti- all topologically different types of neutrino events lead to 

■ cally reviewed and details are in the following: the unified numerical oscillation parameters, say, Am 2 = 

2.4 x lCT 3 eV 2 and sin 2 20 = 1.0 0. 

(1) During 1980's Kamiokande and 1MB observed smaller However, these parameters are obtained from the anal- 
atmospheric neutrino flux ratio Vlt /v e than the pre- ysis of the zenith angle d i stribution f the incident neutri- 
dicted value [lj . nos w hi cn are based on the survival probability of a given 

(2) Kamiokande found anomaly in the zenith angle distri- flavor ( see Eq.(7)). The most important result among the 
bution z . achievements on neutrino oscillation made by Super- 

(3) Super-Kamiokande found v il -v T oscillation [3] and Kamiokande Collaboration is the finding of the maximum 
Soudan2 and MACRO confirmed the Super-Kamiokande oscillation in neutrino oscillation, because it is the ulti- 
result [4J. mate result in the sense that they observe the oscillation 

(4) K2K, the first accelerator-based long baseline experi- pat tern itself directly in neutrino oscillation. 

ment, confirmed atmospheric neutrino oscillation [SI. mi- <■ ■ ■ ■ 

(5) MINOS's precision measurement gives the consistent J^ f cc ° unt of .f u fa f; tors ^T a ^7 e ' :t 1S 

, , . , , c T ^ • i i rn\ natural that the maionty believes the lindmg ot the u — t 

results with Super-Kamiokande ones lol. . ... im * j; T . . . . ° ,. . ,. 

neutrino oscillation by Super- Kamiokande Collaboration. 

It is well known that Super-Kamiokande Collaboration However, it should be emphasized strongly that Super- 
examined all possible types of the neutrino events, such Kamiokande Collaboration put the fundamental assump- 
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tion in all possible analyses of the atmospheric neutrino 
oscillation which is never self-evident and should be care- 
fully examined. This assumption is that the directions of 
the incident neutrinos are approximately the same as those 
of emitted leptons. 

In order to avoid any misunderstanding toward the SK 
assumption on the direction, we reproduce this assump- 
tion from the original SK papers and their related papers 
in italic. 

[A] Kajita and Totsuka [8] state 0: 

"However, the direction of the neutrino must be 
estimated from the reconstructed direction of the 
products of the neutrino interaction. In water Cheren- 
kov detectors, the direction of an observed lepton is 
assumed to be the direction of the neutrino. Fig. 11 
and Fig. 12 show the estimated correlation angle be- 
tween neutrinos and leptons as a function of lepton 
momentum. At energies below 400 MeV/c, the lep- 
ton direction has little correlation with the neutrino 
direction. The correlation angle becomes smaller with 
increasing lepton momentum. Therefore, the zenith 
angle dependence of the flux as a consequence of 
neutrino oscillation is largely washed out below 400 
Me V/c lepton momentum. With increasing momen- 
tum, the effect can be seen more clearly." 

[B] Ishitsuka [5] statefl 

" 8.4 Reconstruction of L v 

Flight length of neutrino is determined from the 
neutrino incident zenith angle, although the energy 
and the flavor are also involved. First, the direction 
of neutrino is estimated for each sample by a dif- 
ferent way. Then, the neutrino flight lenght is cal- 
culated from the zenith angle of the reconstructed 
direction. 

8.4-1 Reconstruction of Neutrino Direction 
FC Single-ring Sample 

The direction of neutrino for FC single-ring sam- 
ple is simply assumed to be the same as the recon- 
structed direction of muon. Zenith angle of neu- 
trino is reconstructed as follows: 



cos 0^ ec — cos &„ 



.17) 



,where cos(9^ ec and cos 0^ are cosine of the recon- 
structed zenith angle of neutrino and muon, respec- 
tively.'" 

[C] Jung, Kajita et al. []Ij] state 0: 

"At neutrino energies of more than a few hundred 
Me V, the direction of the reconstructed lepton ap- 
proximately represents the direction of the original 
neutrino. Hence, for detectors near the surface of 
the Earth, the neutrino flight distance is a function 



see page 101 of the paper concerned. 

2 see page 138 of the paper concerned. 

3 ._„,.., 



sec 



page 453 of the paper concerned. 



of the zenith direction of the lepton. Any effects, 

such as neutrino oscillations, that are a function 

of the neutrino flight distance will be manifest in 

the lepton zenith angle distributions." 

Hereafter, we call the fundamental assumption by the 
Super-Kamiokande Experiment simply as the SK assump- 
tion on the direction. 

Among various types of the neutrino events analyzed 
by Super-Kamiokande, the most important events are sin- 
gle ring muon(electron) events which are generated in the 
detector and terminate in it, because these events give 
more essential physical quantities for clear interpretation 
on neutrino oscillation, namely, the kinds of the neutrino, 
the transferred energies and the directions of the produced 
leptons. These single ring muon events are generated from 
the quasi elastic scattering (QEL). If the existence of neu- 
trino oscillation is verified surely in the analyses of single 
ring muon events among Fully Contained Events under the 
SK assumption on the direction, then we can expect the 
corroboration of the oscillation in the analyses of other 
types of the interaction with less accuracy. Therefore, the 
SK assumption on the direction should be carefully ex- 
amined in the analysis of single ring events due to QEL 
among Fully Contained Events. 

Our paper is organized as follows. In section 2, we treat 
the differential cross section for QEL in the stochastic 
manner as exactly as possible and obtain the zenith an- 
gle distributions of the emitted leptons, cos 9^, for given 
neutrinos with definite zenith angles, taking account of 
the azimuthal angles of the emitted leptons in QEL. As 
a result of it, we show that the SK assumption on the di- 
rection does not hold any more for the incident neutrinos 
with smaller energies. In section 3, we examine the SK as- 
sumption on the direction in the light of L v and L^ and 
reach the same conclusion obtained in the section 2, as it 
must be. In section 4, we examine all possible choices of 
combination of L and E, namely, L v /E v , L v / E^, L^/E v , 
L^/E^, in L/E analysis and show that the maximum os- 
cillation can be realized in the combination of L v j E v , as 
it must be. In section 5, we compare our results obtained 
from the numerical experiments with real observation ob- 
tained by Super-Kamiokande Collaboration. In section 6, 
we conclude that SK cannot observe the maximum oscil- 
lation in their L^/ E v analysis. 

Here, we designate the minimum extrema in L/E dis- 
tribution for neutrino oscillation as the maximum oscil- 
lation by the terminology already utilized in the Super- 
Kamiokande Collaboration [T3"] . 

2 Single Ring Events among Fully Contained 
Events which are Produced by Quasi Elastic 
Scsattering. 

2.1 Differential cross section of quasi elastic 
scattering and influence over various quantities 
concerned 

As stated in Introduction, the finding of the observation 
of the maximum oscillation in the L/E analysis is the ulti- 
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mate verification of the finding of the neutrino oscillation 
by Super-Kamiokande. For the examination of the Super- 
Kamiokande's assertion, we analyze the L/E distribution 
of the single ring events among Fully Contained Events. 

In order to examine the validity of the SK assumption 
on the direction, we consider the single ring events due to 
the following quasi elastic scattering(QEL): 



v e + n 
v^ + n - 
D e +p - 



l'n 



P 



-> p + e 

>P + P 

* n + e 

n + [i 



+ 



(1) 



The differential cross section for QEL is given as follows 

nu. 



d&e(e)( E Hv) 



dQ 2 




where 



A{Q 2 ) = % 



,2/ Q 2 



A I M 2 



ff AQ2 



+*(£.-*. 



\M 2 AM 4 
B(Q 2 ) = (h + / 2 ).9i0 2 , 



9?( 4 



SI 

M 2 



C(Q 2 



M 2 



f 2 , f 2 g 

1 2 4M 2 



■<71 2 



The signs + and — refer to 



V(e) 



and Vftte) for charged 



current (c.c.) interactions, respectively. The Q 2 denotes 
the four momentum transfer between the incident neu- 
trino and the charged lepton. Details of other symbols are 
given in |llj . 

The relation among Q 2 , E v rp\, the energy of the in- 
cident neutrino, Ei, the energy of the emitted charged 
lepton (muon or electron or their anti-particles) and (9 S , 
the scattering angle of the emitted lepton, is given as 

Q 2 = 2E v(B) Et(l - cos6 s ). (3) 

Also, the energy of the emitted lepton is given by 



Ei — E v {p) — 



2M' 



(4) 



Now, let us examine the magnitude of the scattering 
angle of the emitted lepton in a quantitative way, as this 
plays a decisive role in determining the accuracy of the di- 
rection of the incident neutrino, which is directly related 
to the reliability of the zenith angle distribution of sin- 
gle ring muon (electron) events in the Super-Kamiokande 
Experiment. 




E[i/Ev 
Fig. 1. Relation between the energy of the muon and its scat- 
tering angle for different incident muon neutrino energies, 0.5, 
1, 2, 5, 10 and 100 GeV. 



0.050 



0.040 - 



« 0.030 
T3 0.020 



0.010 



Ev = 2 GeV 
I \ 1 GeV 












0.5 GeV 









20 40 60 80 100 120 140 160 IS 

0s ( degree ) 



Fig. 2. Distribution functions for the scattering angle of the 
muon for muon-neutrino with incident energies, 0.5 , 1.0 and 
2 GeV. Each curve is obtained by the Monte Carlo method 
(one million sampling per each curve). 



By using Eqs. (0) to $$, we obtain the distribution 
function for the scattering angle of the emitted leptons 
and the related quantities by a Monte Carlo method. The 
procedure for determining the scattering angle for a given 
energy of the incident neutrino is described in Appendix 
A. Figure Q] shows this relation for muon, from which we 
can easily understand that the scattering angle 8 S of the 
emitted lepton ( muon here ) cannot be neglected. For a 
quantitative examination of the scattering angle, we con- 
struct the distribution function for 9 S of the emitted lepton 
from Eqs. ([2J to (Ql by using the Monte Carlo method. 

Figure [5] gives the distribution function for 6 S of the 
muon produced in the muon neutrino interaction. It can 
be seen that the muons produced from lower energy neu- 
trinos are scattered over wider angles and that a consid- 
erable part of them are scattered even in backward direc- 
tions. Similar results are obtained for anti-muon neutrinos, 
electron neutrinos and anti-electron neutrinos. 

Also, in a similar manner, we obtain not only the dis- 
tribution function for the scattering angle of the charged 
leptons, but also their average values < 9 S > and their 
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standard deviations a s . Table 1 shows them for muon neu- 
trinos, anti-muon neutrinos, electron neutrinos and anti- 
electron neutrinos. From Table 1, it seems to be clear that 
the scattering angles could not be neglected, taking ac- 
count of the fact that the frequency of the neutrino events 
with smaller energies is far larger than that of the neu- 
trino events with larger energies due to high steep of the 
neutrino energy spectrum. However, Super-Kamiokande 
Collaboration assume that the direction of the neutrino is 
approximately the same as that of the emitted lepton even 
for the neutrino events with smaller energies, as cited in 
the three passages mentioned above [5], [S])[ID]- However, 
it has never been verified by Super-Kamiokande Collabo- 
ration. 



2.2 Influence of azimuthal angle in QEL over the 
zenith angle of single ring events 



In the present subsection, we examine the effect of the 
azimuthal angles, 0, of the emitted leptons over their own 
zenith angles, 0^(ji)) , for given zenith angles of the incident 
neutrinos, Q v (p)) m QEL, which was not be considered 
in the detector simulation carried by Super-Kamiokande 
Collaboration^. The influence of this effect over the zenith 
angle cannot be neglected particularly in horizontal-like 
neutrino events. 

For three typical cases (vertical, horizontal and diag- 
onal), Figure [3] gives a schematic representation of the 
relationship between, #„( P ), the zenith angle of the inci- 
dent neutrino, and ( 9 S , </>), a pair of scattering angle of 
the emitted lepton and its azimuthal angle. 

From Figure |3}a, it can been seen that the zenith an- 
gle O^p.) of the emitted lepton is not influenced by its 



in the vertical incidence of the neutrinos 



Mp) 



0°), as 



it must be. From Figure [3}b, however, it is obvious that 
the influence of <fi of the emitted leptons on their own 
zenith angle is the strongest in the case of horizontal in- 
cidence of the neutrino (6 v (p) = 90°). Namely, one half of 
the emitted leptons are recognized as upward going, while 
the other half is classified as downward going ones. The 
diagonal case ( Q v (p) = 43°) is intermediate between the 
vertical and the horizontal. In the following, we examine 
the cases for vertical, horizontal and diagonal incidence 
of the neutrinos with different energies, say, -E„(p) = 0.5 



GeV, E, 
cases. 



u{u) 



1 GeV and E, 



u{9) 



5 GeV, as the typical 



4 Throughout this paper, we measure the zenith angles of 
the emitted leptons from the upward vertical direction of the 
incident neutrino. Consequently, notice that the sign of our 
direction is opposite to that of the Super-Kamiokande Exper- 





a. 6v = 0° 



b. 9v = 90° 




c. 9v = 43° 



Fig. 3. Schematic view of the zenith angles of the charged 
muons for different zenith angles of the incident neutrinos, fo- 
cusing on their azimuthal angles. 



2.3 Dependence of the spreads of the zenith angle for 
the emitted leptons on the energies of the emitted 
leptons for different incident directions of the 
neutrinos with different energies 



The detailed procedure for the Monte Carlo simulation is 
described in Appendix A. We give the scatter plots be- 
tween the fractional energies of the emitted muons and 
their zenith angle for a definite zenith angle of the inci- 
dent neutrino with different energies in Figures [4] to [6J In 
Figure 21 we give the scatter plots for vertically incident 
neutrinos with different energies 0.5,1 and 5 GeV. In this 
case, the relations between the emitted energies of the 
muons and their zenith angles are unique, which comes 
from the definition of the zenith angle of the emitted lep- 
ton. However, the densities (frequencies of the event num- 
ber) along each curves are different in position to position 
and depend on the energies of the incident neutrinos. Gen- 
erally speaking, densities along the curves become greater 



iment ( our cos^,,^) = - cost 



u(i>) 



in SK) 



toward cos# 



fj-(p) 



1. In this case, cosi 



V0*) 



is never in- 
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Table 1. The average values < 9 B > for scattering angle of the emitted charged leptons and their standard deviations o B for 
various primary neutrino energies E v ^y 



E V{D) (GeV) 


angle 
(degree) 


^M(A) 


u v(n) 


v e 


^e 


0.2 


<e s > 


89.86 


67.29 


89.74 


67.47 


o B 


38.63 


36.39 


38.65 


36.45 


0.5 


<9 B > 


72.17 


50.71 


72.12 


50.78 


o B 


37.08 


32.79 


37.08 


32.82 


1 


<9 B > 


48.44 


36.00 


48.42 


36.01 


o B 


32.07 


27.05 


32.06 


27.05 


2 


<9 B > 


25.84 


20.20 


25.84 


20.20 


o B 


21.40 


17.04 


21.40 


17.04 


5 


<9 B > 


8.84 


7.87 


8.84 


7.87 


o B 


8.01 


7.33 


8.01 


7.33 


10 


<9 S > 


4.14 


3.82 


4.14 


3.82 


o B 


3.71 


3.22 


3.71 


3.22 


100 


<9 B > 


0.38 


0.39 


0.38 


0.39 


o B 


0.23 


0.24 


0.23 


0.24 



flucnced by the azimuthal angel in the scattering by the 
definition Q 

On the contrast, it is shown in Figure [S] that the hor- 
izontally incident neutrinos give the widest zenith angle 
distributions for the emitted muons with the same frac- 
tional energies due to the effect of the azimuthal angles. 
The lower the energies of the incident neutrinos are, the 
wider the spreads of the scattering angles of emitted muons 
Ofj, become, which leads to wider zenith angle distribu- 
tions for the emitted muons. As easily understood from 
Figure |51 the diagonally incident neutrinos give the inter- 
mediate zenith angle distributions for the emitted muons 
between those for vertically incident neutrinos and those 
for horizontally incident neutrinos. 

In Figures [7] to [9j we express Figures |4] to \6\ in a differ- 
ent way. We sum up muon events with different emitted 
energies for given zenith angles. As the result of it, we 
obtain frequency distribution of the neutrino events as a 
function of cosO^ for different incident directions and dif- 
ferent incident energies of neutrinos. 

In Figures [TJa) to [7Jc) , we give the zenith angle dis- 
tributions of the emitted muons for the case of vertically 
incident neutrinos with different energies, say, E v = 0.5, 
1 and 5 GeV. 

Comparing the case for 0.5 GeV with that for 5 GeV, 
we understand the big contrast between them as for the 
zenith angle distribution. The scattering angle of the emit- 
ted muon for 5 GeV neutrino is relatively small (See, Ta- 
ble 1), so that the emitted muons keep roughly the same 
direction as their original neutrinos. In this case, the ef- 
fect of their azimuthal angle on the zenith angle is also 
smaller. However, in the case for 0.5 GeV which is the 
dominant energy for single ring muon events in the Super- 



5 The zenith angles of the particles concerned are measured 
from the vertical direction. 



Kamiokande, there is even a possibility for the emitted 
muon to be emitted in the backward direction due to the 
larger angle scattering, the effect of which is enhanced by 
their azimuthal angle. 

The most frequent occurrence in the backward direc- 
tion of the emitted muon appears in the horizontally inci- 
dent neutrino as shown in Figs. 8(a) to 8(c). In this case, 
the zenith angle distribution of the emitted muon should 
be symmetrical with regard to the horizontal direction. 
Comparing the case for 5 GeV with those both for ^0.5 
GeV and for ~1 GeV, even 1 GeV incident neutrinos lose 
almost the original sense of the incidence if we measure 
it by the zenith angle of the emitted muon. Figures IHJa) 
to [9jc) for the diagonally incident neutrinos tell us that 
the situation for diagonal case lies between the case for the 
vertically incident neutrinos and that for horizontally inci- 
dent ones. SK in the figures denotes the SK assumption on 
the direction of incident neutrinos. From the Figures EJa) 
to injc) , it is clear that the scattering angles of emitted 
muons influence their zenith angles, which is enhanced by 
their azimuthal angles, particularly for more inclined di- 
rections of the incident neutrinos. 



3 Super-Kamiokande Assumption on the 
Direction in the Light of L u and L M 

In the previous section, we show that the SK assumption 
on the direction does not hold as for scattering angles of 
the leptons even if statistically. This assumption is logi- 
cally equivalent to the statement that L v is approximately 
the same as L M in L/E analysis , where L v denotes the dis- 
tance on the incident neutrino from the interaction point 
of the neutrino events to the intersection of the Earth sur- 
face toward its arriving direction and L^ denotes the cor- 
responding distance on the emitted muon. Consequently, 
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Fig. 4. The scatter plots between the fractional energies of the produced muons and their zenith angles for vertically incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Fig. 5. The scatter plots between the fractional energies of the produced muons and their zenith angles for horizontally incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Fig. 6. The scatter plots between the fractional energies of the produced muons and their zenith angles for diagonally incident 
muon neutrinos with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Fig. 7. Zenith angle distribution of the muon for the vertically incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stands for the corresponding ones under the SK assumption. 
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Fig. 8. Zenith angle distribution of the muon for the horizontally incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stands for the corresponding ones under the SK assumption. 
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Fig. 9. Zenith angle distribution of the muon for the diagonally incident muon neutrino with 0.5 GeV, 1 GeV and 5 GeV, 
respectively. The sampling number is 10000 for each case. SK stands for the corresponding ones under the SK assumption. 
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Fig. 10. Schematic view of relations among L„, L M , 6 S and 
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In section 3.1 we give the procedure how to obtain 
L^ from a neutrino event with given L v in the stochastic 
manner. In section 3.2 we give the correlations between 
L v and L^, taking account of the effect of the backscat- 
tering as well as the effect of the azimuthal angle in the 
QEL in stochastic manner. As the result of it, we show 
that L v w L^, namely the SK assumption on the direc- 
tion, does not hold even if statistically in both the ab- 
sence and the presence of neutrino oscillation (Figure [T2l 
and Figure H"3|) . Also, we treat the correlation between 
E v and E^, in the stochastic manner. We show that the 
approximation of E v with E^ by Super-Kamiokande Col- 
laboration does not make so serious error compared with 
the approximation of L v by L^, although their treatment 
is theoretically unsuitable (Figure [5]). 

In section 4, we show that L v jE v distribution can 
reproduce the minimum extrema for neutrino oscillation 
which SK's neutrino oscillation parameters demand and 
, furthermore, it may give the differnt mimimum extrema 
in the neutrino oscillation under the different neutrino 
oscillation parameters from SK's. We show LyjE v dis- 
tribution can reproduce the minimum extrema for neu- 
trino oscillation (the maximum oscillation) which Super- 
Kamiokande Collaboration demand, by using their neu- 
trino oscillation parameters (Am 2 = 2.4 x 10 -3 eV 2 and 



2 2H 



1.0). Furthermore, it may give the different 



minimum extrema in the neutrino oscillation under the 
different neutrino oscillation parameters from the Super- 
Kamiokande Collaboration. This fact denotes that our nu- 
merical computer experiment is done in a correct manner. 



<p s sampling 
determining Z„ 



sampling of a 

interaction point 

judging 

FC or PC 



(L V L U ) for FC 



GEANT3 



Fig. 11. The procedure for our numerical experiment for ob- 
taining L M from a given L v . 



if our indication on the invalidity of the SK assumption on 
the direction is correct, the same conclusion should be ex- 
pected in the relation between L v and L M . In the present 
section and subsequent sections, we examine directly the 
validity of the implicit SK assumption that L v is approx- 
imated by L^, taking into consideration the neutrino en- 
ergy spectrum at the Super-Kamiokande site. 



3.1 Derivation of L^ from a given L v in the single ring 
muon event among Fully Contained Events 

In our numerical computer experiment, we obtain single 
ring muon events among the Fully Contained Events re- 
sulting from QEL in the virtual SK detector, the details 
of which are described in Appendix A. For the neutrino 
event with a definite neutrino energy thus generated, we 
simulate its interaction point inside the detector and the 
emitted energy of the muon concerned which gives its scat- 
tering angle uniquely. The determination of the neutrino 
energy, the emitted energy of the muon and its scattering 
angle are described in Appendix A. The muon thus gener- 
ated is pursued in the stochastic manner by using GEANT 
3 and finally we judge whether the muon concerned stops 
inside the detector (the Fully Contained Event) or passes 
through the detector (the Partially Contained Event). For 
Fully Contained Events thus obtained, we know the direc- 
tions of the incident neutrinos, the generation points and 
termination points of the events generated inside the de- 
tector, the emitted muon energies, their scattering angles 
and their azimuthal angles in QEL which give their zenith 
angles, L„ and L M finally. 

In Figure ITUl we show the relation between L v and L^ 
schematically. Figure [TT] shows the procedure for obtain- 



6 The azimuthal angle is but that in QEL, not that with 
regard to the Earth here. 
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Fig. 12. Correlation diagram for L v and L M without oscillation 
for 1489.2 live days. The blue points and orange points denote 
neutrino events and ani-neutrino events, respectively. 
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Fig. 13. Correlation diagram for L„ and L M with oscillation 
for 1489.2 live days. The blue points and orange points denote 
neutrino events and ani-neutrino events, respectively. 



ing L/j, from L v which is equivalent to the corresponding 
procedure for obtaining cos9^ from cos9 v . 

The relation between direction cosine of the incident 
neutrino, (4(5)i"ii/(5))n v (5)), and that of the correspond- 
ing emitted lepton, (£ r , m v , n r ) , for a given scattering an- 
gle, a , and its azimuthal angle, 4>, resulting from QEL is 
given in Appendix A. 

L v and L^ are functions of the direction cosine of the 
incident neutrino, cos9 v , and that of emitted muon, cos9^, 
respectively and they are given as, 



L v =R g x (r SK cos9„ + \ r 2 SK cos 2 9 v + 1 - r 2 SK ) (5 - 1) 



Rg x (rsKCOs9 fJ 



SK 



COS' 



' SK) 



(5-2) 



Ev-Eu correlation for FC events 

1489.2 days Oscillation Am 2 = 2.4* 10~ 3 eV 2 




2 

Ev ( GeV ) 

Fig. 14. The correlation diagram between E v and _E M for os- 
cillation for 1489.2 live days. The continous line denotes the 
polynomial expression by Super-Kamiokande Collabolation. 



where R g is the radius of the Earth and r$K = 1 — 
Dsk/Rq, with the depth, D$k, of the Super-Kamiokande 
Experiment detector from the surface of the Earth. It 
should be noticed that the L v and L M are regulated by 
both the energy spectrum of the incident neutrino and 
the production spectrum of the muon (QEL in the present 
case) . Consequently, their mutual relation is influenced by 
either the absence of the oscillation or the presence of the 
oscillation which depend on the combination of the oscil- 
lation parameters. 



3.2 The correlation between L v and L M 

In Figure 1121 we give the correlation diagram between L v 
and L^ for single ring muon events among Fully Contained 
Events for the 1489.2 live days in the absence of neu- 
trino oscillation which corresponds to the actual Super- 
Kamiokande Experiment [7. In Figure [T2l blue points de- 
note neutrino events while orange points denote anti- 
neutrino events. Throughout all correlation diagrams in 
the present paper, blue points and orange ones have the 
same meaning shown in Figure [T2] The aggregates of the 
(anti-) neutrino events which correspond to a definite com- 
bination of L v and L^ are essentially classified into four 
groups in the following: 

Group A is defined as the aggregate for neutrino events 
in which both L v and L M are rather small. It denotes that 
the downward neutrinos produce the downward muons 
with smaller scattering angles. In this case, the energies of 
the produced muons are near the energies of the incident 
neutrinos due to smaller scattering angles. 

Group B is defined as the aggregate for neutrino events 
in which both L v and L^ are rather large. It denotes 
that the upward neutrinos produce upward muons with 
smaller scattering angles. In this case, the energy relation 
between the incident neutrinos and the produced muons 
is essentially the same as in Group A, because the flux 
of the upward neutrino events is symmetrical to that of 
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the downward neutrino events in the absence of neutrino 
oscillation. 

Group C is defined as the aggregate for neutrino events 
in which L v are rather small and L M are rather large. 
It denotes that the downward neutrinos produce the up- 
ward muons by the possible effect reusulting from both 
backscattering and azimuthal angle in QEL. In this case, 
the energies of the produced muons are smaller than those 
of the energies of the incident neutrinos due to larger scat- 
tering angles. 

Group D is defined as the aggregate for the neutrino 
events in which L v are rather large and L^ are rather 
small. It denotes that the upward neutrinos produce the 
downward muons. The energy relation between the inci- 
dent neutrinos and the produced muons is essentially the 
same as in Group C in the absence of neutrino oscillation. 

It is clear from Figure [T^] that there exist the symme- 
tries between Group A and Group B, and also between 
Group C and Group D, which reflect the symmetry be- 
tween the upward neutrino flux and the downward neu- 
trino one for null oscillation. 

In Figure 1131 we give the correlation between L v and 
L^ under their neutrino oscillation parameters, say, Am 2 — 
2.4 x 10~ 3 eV 2 and sin 2 29 = 1.0 7 . In the presence of 
neutrino oscillation, the property of the symmetry which 
holds in the absence of neutrino oscillation (see (Group 
A and Group B) and/or (Group C and Group D) in Fig- 
ure [T2]) is lost due to the different incident neutrino fluxes 
in the upward direction and downward one. If we com- 
pare Group A with Group B, the event number in Group 
B (upward v — > upward /i) is smaller than that in group A 
(downward v — > downward /i), which comes from smaller 
flux of the upward neutrinos. The similar relation between 
Group C (downward v — > upward /i ) and Group D (up- 
ward v —¥ downward \x) is held in Figure 1131 

Summarizing the characteristics among Groups A to D 
in the Figures IT21 and ITUl we could conclude that (Group 
A and Group B) and (Group C and Group D) are in sym- 
metrical situations in the absence of neutrino oscillation, 
while such a symmetrical situation is lost in the presence 
of neutrino oscillation. Also, it is clear from Figures H"2l 
and 1131 that L v sa L M , namely the SK assumption on the 
direction, does not hold both in the absence of neutrino 
oscillation and in the presence of neutrino oscillation even 
if statistically. 

Here, it should be noticed that the approximation of 
L v k, L^ does not hold completely in the region C and 
region D. The event numbers in Group C and Group D 
could not be neglected among the total event number con- 
cerned. In these regions, neutrino events consist of those 
with backscattering and/or neutrino events in which the 
neutrino directions are horizontally downward (upward), 
but their produced muons turn upward (downward) re- 
sulting from the effect of azimuthal angles in QEL. 



simulation, by the following equation[5] (see page 135 of 
the paper concerned) : 



E, 



v,SK 



E^ x (a + b x x + c x x 2 + d x a; 3 ), (6) 



where x = logioiE^). 

The idea that E v could be approximated as the poly- 
nomial means that there is unique relation between E v 
and E^. However, in the light of stochastic characters in- 
herent in both the incident neutrino energy spectrum and 
the production spectrum of the muon, such a treatment 
is not suitable theoretically, which may kill a real correla- 
tion effect between the incident neutrino energy and the 
emitted muon energy. In Figure Q31 we give the correlation 
between E v and E^ together with that obtained from the 
polynomial expression by Super-Kamiokande Collabora- 
tion under their neutrino oscillation parameters and their 
incident neutrino energy spectrum|12 . It is clear from the 
figure that the part of the lower energy incident neutrino 
deviates largely from the approximated formula, which re- 
flects explicitly the stochastic character of QEL. We can 
give the similar relation for null oscillation, the shape of 
which may be different from that with oscillation due to 
the difference in the incident neutrino energy spectrum. 

Thus, we could choose four combinations, namely L v / E v , 
L^/ E^, L^/ E v and L v /E^ for the examination of maxi- 
mum oscillations due to neutrino oscillation in L/E anal- 
ysis. However, only the combination of L^/E^ out of these 
four combinations can be physically measurable. 



4 Summary 

Since one cannot measure L v and E v , so one is forced to 
utilize L M and E^ in the L/E analysis in place of them. 
Then, Super-Kamiokande Collaboration assume that the 
direction of the incident neutrino is the same as that of 
the emitted lepton the SK assumption on the direction and 
E v can be estimated from the some polynomial formula of 
the variable E^ in L/E analysis. However, it is clear from 
Figures 12 and 13 that the SK assumption on the direction 
does not hold even approximately and the transformation 
of E^ into E v is not uniquely. 

In the Part 2 of the subsequent paper, we apply the 
results from Figures 12, 13 and 14 to L/E analysis and 
conclude that one cannot obtain the maximum oscillation 
in L/E analysis which shows strongly the oscillation pat- 
tern from the neutrino oscillation. 



3.3 The correlation between E v and E u 



Super-Kamiokande Collaboration estimate E v from E^, 
the visible energy of the muon, from their Monte Carlo 
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A Monte Carlo Procedure for the Decision 
of Emitted Energies of the Leptons and Their 
Direction Cosines 

Here, we give the Monte Carlo Simulation procedure for 
obtaining the energy and its direction cosines, (l r , m r , n r ), 
of the emitted lepton in QEL for a given energy and its 
direction cosines, (l,m,n), of the incident neutrino. 

The relation among Q 2 , E u t 9 ), the energy of the in- 
cident neutrino, -EW), the energy of the emitted lepton 
(muon or electron or their anti-particles) and 6 S , the scat- 
tering angle of the emitted lepton, is given as 



Q 2 = 2E v{D) E m (l -cos9 s ). 
Also, the energy of the emitted lepton is given by 



E i(i) — Ev{v) - 



91 

2M' 



(A-l) 



(A-2) 



Procedure 1 

We decide Q 2 from the probability function for the differ- 
ential cross section with a given E v i&\ (Eq. © in the text) 
by using the uniform random number, £, between (0,1) in 
the following 



cf 



Pe(e)(K { ,),Q 2 )dQ 2 



(A-3) 



where 



P t(J){ E v(u),Q ) - 

dQ^~ 



(E^),Q 2 ) 




dQ 2 . 



(A-4) 



From Eq. (A-l), we obtain Q 2 in histograms together 
with the corresponding theoretical curve in Figure[T5l The 
agreement between the sampling data and the theoretical 
curve is excellent, which shows the validity of the utlized 
procedure in Eq. (A-3) is right. 

Procedure 2 

We obtain Eug> from Eq. (A-2) for the given E v (p-\ and 
Q 2 thus decided in the Procedure 1. 

Procedure 3 

We obtain cos# s , cosine of the the scattering angle of the 
emitted lepton, for E^g\ thus decided in the Procedure 2 
from Eq. (A-l) . 

Procedure 4 

We decide (j>, the azimuthal angle of the scattering lepton, 




0.4 0.6 
Q 2 /Q 



max 



Fig. 15. The reappearance of the probability function for QEL 
cross section. Histograms are sampling results, while the curves 
concerned are theoretical ones for given incident energies. 




(/, m, n) 



Fig. 16. The relation between the direction cosine of the 
incident neutrino and that of the emitted charged lepton. 
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which is obtained from 



= 2tt£. 



(A-5) 



Here, £ is a uniform random number between (0, 1). 
As explained schematically in the text (see Figure [3] in the 
text), we must take account of the effect due to the az- 
imuthal angle <f> in the QEL to obtain the zenith angle 
distribution both for Fully Contained Events and Partially 
Contained Events correctly. 

Procedure 5 

The relation between direction cosines of the incident neu- 
trinos, {i v (p), rn v ^,n v ^), and those of the corresponding 
emitted lepton, (£ r , rn r , n r ), for a certain 9 S and <f> is given 
as 




/ 



In 



m 



\j£ 2 + m 2 
mn 

yj £ 2 + TTT 2 
V -V£ 2 + m 



Vp 



£v(v\ \ 



,'■' "' ' 1 /sin# s co< 



VI 2 



V(p) 



where n 



v(y) 



cos9 



v{v) 



sin(9 s sir 

\ COS# s 

(A-6) 

and n r = cos#^. Here, 0£ is the 



T-u(p) 



«/ 



zenith angle of the emitted lepton. 

The Monte Carlo procedure for the determination of 9t 
of the emitted lepton for the parent (anti-)neutrino with 
given Q v [p) and £v(p) involves the following steps: 
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We obtain (£ r ,m r ,n r ) by using Eq. (|A-6|) . The n r is 
the cosine of the zenith angle of the emitted lepton which 
should be contrasted to n v , that of the incident neutrino. 
Repeating the procedures 1 to 5 just mentioned above, we 
obtain the zenith angle distribution of the emitted leptons 
for a given zenth angle of the incident neutrino with a def- 
inite energy. 



In the SK analysis, instead of Eq. (|A-6[) . they assume 
n r = n„( P ) uniquely for -E^(p) > 400 MeV. 
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